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ABSTRACT 


This paper reports an attempt to simplify for ease in field use 
the complex numerical model for sound propagation loss prediction 
developed by Dr, C, S, Clay, Captain P, M, Wolff, and Dr, P, R, Tatro 
(LCDR, USN), Due to time required for communications and computa- 
tion, the present results of this model are not available for imme- 
diate use by the operator in the field, 

This computerized nose propagation loss for 
temperature profiles with sea surface temperature varying between 10C 
and 25C and variation in mixed layer depth from the surface to 300 
meters, From these data formulas and graphs were developed to predict 
convergence zone range, width, and intensity and the 100—decibel propa- 
gation loss range, These results are intended for operational use by 


field personnel, 
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CHAPTR I 


INTRODUCTION 


A numerical model for calculating sound propagation loss from a 
source at a given depth to a receiver at a given depth in the sea is 
presently employed at the Fleet Numerical Weather Facility, (FNWF), 
Monterey, California, This model is used to provide tactical infor- 
mation to operational naval forces and for other strategic uses, How- 
ever, due to turnaround time required for computation and communica- 
tions, the results are not available for immediate operational usage 
in the “field, 

This model was developed and adapted for computers in November 
1967 by C, S, Clay, P, M, Wolff, and P, R, Tatro, Recent comparisons 
of predicted values with values obtained in controlled field studies 
have shown the model to produce excellent results, The Bolt, Beranek, 
and Newman (1966) report gave finer detail on sound loss, but the FNWF 
model compared very favorably over the entire sound propagation loss 
path, 

Having worked operationally in this field as ASW officer on the 
USS VESOLE (DDR-878) and on the USS RICHARD E, BYRD (DDG-23), the author 
has firsthand experience in the many frustrations encountered in 
coping with sound detection in the sea, This frustration is shared by 
other operators as noted by Glennon (1965), By using the results of 
the FNWF model as a standard, it is intended to develope some simple 
equations and graphical displays for use by the operator directly in 
the field, 

This work is concentrated on the convergence zone range, width, 


and intensity for convergence zones one, two, and three, Data and graphs 


9 


are also presented for the 100-decibel propagation loss range, Addi- 
tional runs would be required to utilize the 100-decibel data in a more 
meaningful manner, 

This paper is divided into three principal parts, The first 
part is an explanation of the basic features of the sound propagation 
loss model used to develop the data, Since the model description has 
not been published, this information is derived from personal consulta- 
tion with Dr, P, R, Tatro and Dr, C, S, Clay, 

The second part will include an approach to obtaining data for 
reduction, and a description of the methods used for reducing the 
data obtained, The fine atlas of Pacific water conditions developed 
by Muromtsev (1963) has been of considerable help in developing this 
ocean model, 

The last section will be the display of graphical results and the 
equations developed for operational use, These calculations will be 
shown with the realization that some differences may exist between 
these results and actual conditions in the field, It is intended that 
the model predictions will be accurate enough to be of some operational 


use 
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CHAPTER II 


SOUND PROPAGATION LOSS MODEL 


The sound propagation loss model was developed to calculate the 
transmission loss for low frequency sound in the world oceans, Losses 
developed in the model depend on sound velocity profile, depth, sea 
state, sound frequency, and the type of ocean bottom present, These 
parameters in turn are dependent on such variables as sea surface 
temperature, mixed layer depth, temperature profile, source and re- 
ceiver depths, bottom roughness, and salinity, The model was developed 
in such a manner that these variables may be changed independently as 
demanded by environmental conditions or operational need, Anderson 
and Lesser (1959) indicate some effects of these changes, 

The procedures in the construction of this model have been used for 
experimental analysis by various research facilities, The model com- 
bines these procedures with modern high speed computers to provide 
products for sonar systems on a worldwide scale, 

Sound loss in the model is assumed to occur over ten different 
types of paths in the ocean, These paths are; 

(1) First bottom bounce path up 

(2) First bottom bounce path down 

(3) Second bottom bounce path up 

(4) Second bottom bounce path down 

(5) Convergence zone path up 

(6) Convergence zone path down 

(7) Direct path up 

(8) Direct path down 

(9) Surface duct path 


ie 


(10) Subsurface duct path, 

Urick (1964) and Anderson and Pedersen (1965) have further amp- 
lifying details on these transmission paths in the ocean, 

For further clarification of these paths, a graphical display of 
each path is shown in Figures 1, 2, 3, and 4, 

Sound energy 1s traced along ray paths for every two degrees 
from horizontal (0 degrees) to nearly vertical (88 degrees), Sound 
propagation loss is computed and stored in the computer at 1000 yard 
(0,5 nautical mile) increments for ranges extending from 1000 yards 
to 250,000 yards, Since these are one-way transmission losses, they 
should be doubled if used with active sonar systems, 

Sound energy loss contributions to the total energy loss along 
the various paths are calculated s crcageie ag the total loss is the 
ratio of total energy available via the various paths to an initial 
energy assumed to be unity, This is conventionally expressed as a 
logarithmic unit, the decibel, These components are independent of 
each other and the method for obtaining each component will be covered 


in the following paragraphs, 


Losses Along Direct Path, Convergence Zone and Bottom Reflected Paths 


Losses via direct path, convergence zone, and bottom bounce are 


calculated by use of the following equation; 


= ss 
Loss = 10 log (——< 5;) + 10 log LetheaRe+ Byer ily (1) 
where 
R = range from source to receiver 
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Figure 4 
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L = vertical distance between spreading rays in the range- 
depth plane 
a = absorption coefficient 
Np = loss due to reflection from bottom 
N, = loss due to reflection from ocean surface 


10 log (——) = cylindrical spreading loss 
cos Oj 


10 log L = differential spreading loss in the range-depth plane, 


Since most of the ocean bottom is considered rough for low fre- 
quency sound, this roughness has been accounted for in the model by 
the term N, in equation (1), Using a varied scale of roughness from 
one to five, where one is the smoothest type bottom and five the 
roughest, tne loss due to bottom variability is computed, The cal- 
culation of the roughness has been made using seismic date, physio- 
graphic provinces, and bathymetric data, 

The composition of the bottom is included in the model through 
N, by using a given propagation loss for certain materials, At present 
only two designations are used, and these factors differentiate between 
the "fast" or low loss bottom of the Atlantic Ocean and the slow’! or 
higher loss bottom encountered in the Pacific Ocean, N, is also a 
function of the angle which the sound ray makes witn the ocean bottom, 

The loss due to surface reflection, N,, is a function of sea 
state, the angle the ray makes with the ocean surface, and the fre- 
quency of the sound traveling through the water, The cylindrical 
spreading loss and the differential spreading loss combine to give the 
geometric spreading loss of the signal, 

Bottom reflected energy is calculated using equation (1) out 
through only two bottom bounces, Energy loss after the second bounce 


is extrapolated using; 
af 





LOSS = Constant + 10 log R+ aR + ABR (2) 
Constant = L-= ( 10 log R, + aR, ) 


L, 


sound energy at end of second bottom bounce 


Rj = range to end of second bottom bounce zone 


R = total range 


eee 


Do 


Lo.) = energy loss at 0 = 0,4 radians 
ay = angle incident ray makes with bottom 


D. = range where no more rays exist that have hit the bottom 
only once, or twice for first or second bottom bounce respec— 


tively, 


Surface Duct and Subsurface Duct Losses 
Surface duct and subsurface duct losses are treated using an 
approximation to a normal mode solution, The equation used for this 


determination is 


—BY =o R 
i = Nag e O e7ah 
A RZo (3) 


Loss = -l10O log i 





Io 
where 
R = range from source to receiver 
I, = initial intensity of signal 


49 = angle of energy trapped in a duct of depth Z 
Xo = half cycle distance 
Y = wavelength of sound 
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a = absorption coefficient 


NS 
tI 


the vertical extent of either the surface or the sub= 
surface duct 


8 and N = empirical constants, 


8 and N may be altered as verification data from field studies 
become available, which permits the model to be "tuned" without diffi- 


culty, 


ay 





CHAPTER I1I 


FI#LD PREDICTION APPROACH 


With this model available for sound loss prediction, the problem 
was to determine the variables to be used to produce information satis-~ 
factory for field prediction, James (1966) and Urick (1967) offer 
guidance in determining the best group of variables, 

The parameters available included sea surface temperature, mixed 
layer depth, temperature gradient, sea state, source depth, receiver 
depth, frequency, bottom depth, bottom type and roughness, and salinity, 
Bach parameter may be varied independently to determine the effects it 
has on the total propagation loss in sound transmission through the 
sea, which is indicated in matrix form in Figure 6, An interpretation 
of the matrix is given in Figure 5, 

The sample output in Figure 6 gives the propagation loss from the 
source to the receiver for a 4000—meter Cige23 feet) water depth, zero 
mixed layer depth, and 15C sea surface temperature, The first, second, 
and third convergence zones are underlined and can be clearly seen in 
the particular example, 

Sea surface temperature, mixed layer depth, and temperature 
gradient have a strong influence on sound propagation and are easily 
determined operationally at the scene; therefore, these variables were 
selected for extensive use in this study, The selection of particular 
ranges of values to be used will be discussed in detail below, 

Primary interest was centered on convergence zones which are 
caused by the partial focusing of refracted ray paths and generally 
require a deep sound channel and large ocean depth, Thus, variation 


in bottom depth over a limited range was included, The focusing 
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Figure 6 


effect can give an intensity from 10 to 10° times the intensity witnout 
this focusing effect, 

Source deptns, receiver depths, and sound frequencies were deter- 
mined through operational considerations, Source depths of 20 feet, 
60 feet, and mixed layer depth plus 100 feet; and receiver depths of 
20 feet, 90 feet, and 300 feet were used, Frequencies used were 
mGO Hz, 3,5 kHz, and 10,0 kHz. 

Sea state was evaluated for states two, four, and six, while 
salinity was assumed constant at 34,70%o. This value was chosen since 
it is the average surface salinity for the Pacific Ocean according to 
Muromtsev (1963). 

The study was confined to the Pacific Ocean and a "slow! or ‘Pil 

bottom type was used, The bottom roughness coefficient used was four, 
since that type bottom occurs most frequently in the Pacific Ocean, 
A bottom depth of 4000 meters was chosen since this approximates the 
depth of much of the water close to the western coast of the United 
States, Data were also obtained for depths of 5000 meters and 6000 
meters; however, these data were not as extensive as those obtained 
mor 1000 meters, 

Mixed layer depths were varied between the surface and 300 meters 
in 100-meter increments, Sea surface temperature was varied between 
10C and 25C, This range of values was chosen as representative of 
most of the Pacific Ocean, Data from Muromtsev (1963) and from per- 
sonnel working extensively with the Pacific Ocean at the FNWF were 
used to determine the initial temperature profiles, Since most of 
the extreme variation in water temperature is observed to occur in the 


upper 400 meters of the ocean, this layer was used in the model for 


2 





strong temperature gradient variations, Sea surface temperatures were 
divided into three groups or families, These groups were 10 to 15C, 
15 to zee, and 20 to, 25G. 

To establish the initial temperature profile for input to the 
sound propagation loss model, the three families were assigned constant 
profiles from 400 meters to the maximum depth of the water (4000 meters, 
5000 meters, or 6000 meters), Examples of each of the three families 
with zero mixed layer depth are shown in Figures 7 and 8, which give 
conditions in the upper 1000 meters, and in the entire 4000 meters, 
respectively, These three families represent broad generalizations of 
the actual conditions in much of the Pacific Ocean, necessary because 
of the scope of this study, 

Using these three families with three temperatures from each 
(tame ly 1s ioyec, 12°5C, 15.0C;. familly 2: 15,00, 17,5C, 20,00; familias 
20,0C, 22,5C, 25,0C) and four mixed layer depths from 0 to 300 meters 
in l0OQ-meter increments, thirty-six initial programs were constructed, 
They were to be evaluated by the computer model using sea states two, 
four, and six; source depths of 20 feet, 60 feet, and mixed layer 
depth + 100 feet; receiver depths of 20 feet, 90 feet, and 300 feet; 
and frequencies of 10,0 kHz, 3,5 kHz, and 100 Hz, 

This initial run produced eighty-one pages of data (sample shown 
on page 22) for each of the thirty-six initial inputs, Upon review of 
the output, it became obvious that meaningful analysis of 250 data 
values for 2,916 data pages would not be possible in the time available, 

Operational usefulness was the main requirement in determining the 
limits of this study, The decision was made to pursue the 100 Hz 


propagation loss profile only, The source depths were reduced to 
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60 feet and mixed layer depth + 100 feet, and the receiver depths 
were limited to 90 feet and 300 fect, It was decided that sea state 
four would give adequate information for operational use, The output 
for each program was thus reduced to four sets of 250 data points, 

More thorough analysis of the initial data showed the necessity 
for additional refinement; it became obvious that more variation would 
be required in sea surface temperature values to obtain meaningful re- 
sults, Therefore, intermediate sea surface temperatures and their 
associated profiles for 11C, 12C, 13C, 14C, 16C, 17C, 18C, and 190 were 
entered into the propagation loss model, No convergence zones appeared 
for sea surface temperatures greater than about 18C at 4000-meter 
water depth and so additional computer runs were made for depths of 
5000 meters and 6000 meters for temperatures from 20 to 25(C, 

Data were developed for the 100-Hz receiver at 90 feet depth with 
a source at 60 feet and at sea state four, These data were examined 
closely for information about convergence zone range, width, and inten- 
sity for convergence zones one, two, and three, The range to the 100- 


decibel propagation loss was also examined, 
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CHAPTER IV 


PROPAGATION LOSS FOR OPERATIONAL USAGE 


The formulas and graphical representations of the data were 
developed to give the clearest and most concise display for opera- 
tional employment, They can be used easily with operator input of sea 


surface temperature, mixed layer depth, and water depth, 


Convergence Zone Range 


The first parameter to be examined was the range to the conver= 
gence zones, Preliminary review indicated the most promising avenue 
of approach related the convergence zone range to sea surface tempera- 
ture and mixed layer depth, Water depth was varied over a sufficient 
range to achieve convergence zone development for sea surface tempera- 
tees up to 25C., 

Data were developed for mixed layer depths of O to 300 meters at 
lOO-meter increments for the temperature limits of 10C and 25C, Due to 
computer time required for each run and data reduction time, effort was 
concentrated on the 100=-meter mixed layer depth for temperatures except 
10 to 15C where computations were carried out for all four mixed layer 
depths and one-degree temperature increments, 

Using the data obtained, graphs were plotted covering the tempera- 
ture band from 10 to 25C for mixed layer depth of 100 meters and for 
all four mixed layer depths from 10 to 15C, A simple equation for 
convergence zone range prediction was developed from these data, The 
development method is shown in the Appendix, This equation uses the 
direct relationship of convergence zone range to sea surface tempera- 
ture and mixed layer depth to predict range to the first, second, and 


third convergence zones, 
28 


R= [ 22,70+0,55(T- 10)] 2 
ee Than 9 0,25zZ1-3¢ a 0.8 0, | M =- 100 
[( ) ( ) ~ (2 ~ 1) (0,5) | Be -A00 (4) 
where 
# 
R = range to the start of the convergence zone (nautical miles) 
T = sea surface temperature (degrees C) 
Z = convergence zone number 


M = mixed layer depth (meters), 


l 


The errors produced by the formula are very small compared with 
the width of the convergence zone (see Tables 4, 5, and 6 in Appendix), 
The average standard deviation of this equation for convergence zones 
one, two, and three at mixed layer depth of 100 meters is 0,61 nautical 
miles, The variance and standard deviation for each of the three con- 


vergence zones is given in Table 1 


RANGE PREDICTION ACCURACY 





Variance (NM) Standard deviation (NM) 
ZONE 1 0,608 0.786 
ZONE 2 OR STN 0.483 
ZONE 3 ora, 0,560 
TABLE I 


Convergence Zone Width 


Width of the convergence zones did not lend itself to any easily- 
predicted formula or organized tabulation of meaningful results when 


compared with mixed layer depth or sea surface temperature, The 
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greatest pattern appeared to be within each group of convergence 

zones, As shown in Table 3, the width of the second convergence zone 
was usually about twice the first and the width of the third convergence 
zone was about three times the first, The width of the first conver- 
gence zone seemed generally to vary between two and four nautical miles; 
a few exceptions were as low as one nautical mile and as high as six 
nautical miles, The second and third convergence zones followed the 
pattern mentioned earlier and were between four and eight nautical miles 
and six and twelve nautical miles respectively, The percentage dis- 
tribution of these data is given in Table 2, for thirty-six observa- 


tions of each convergence zone, 


CONVERGENCE ZONE WIDTH 





ZONK 1 WIDTH (NM) MINIMUM MAXIMUM 
2,0-4,0 42,0 >4.0 an) 6.0 
640% 10% 220% 


ZONE 2 WIDTH (IM) 


4 .0-8,0 </,,0 >8.0 nas 1085 
555% 28 ,0% 16.5% 


ZONE 3 WIDTH (NM) 


6,0-12,0 26,0 >12,0 2,0 14,0 
58.5% 30, 5% 11,0% 
TABLE 2 


This distribution should enable the operator to make a reasonable 


estimate of the zone width, 


Convergence Zone Intensity 


Convergence zone intensity did not lend inself to direct correla- 


tion between mixed layer depth, sea surface temperature, or water depth, 
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of 


Using Figure 9 the operator in the field can immediately assess his 
opportunity of obtaining a target with a source level less than 85 
decibels, between 8) decibels and 95 decibels, and greater than 95 
decibels for each of the first three convergence zones, The break 
points between various decibel levels were determined for operational 
usefulness and have no other distinct meaning, 

The fact that no easily-formulated relationships could be deter- 
mined from the data obtained for convergence —— width and intensity 
does not reduce the importance of the work, It indicates the apparent 
complexity of the convergence zone width and intensity as related to 
sea surface temperature, mixed layer depth, and water depth, However, 
the results obtained are still extremely useful to indicate probability 


of gaining contact by the operator in the field, 


100—Decibel Propagation Loss Level 


The last parameter covered is the range to the 100=decibel propaga=- 
tion loss level, This range was chosen as being representative of the 
capability of a low-frequency passive system, 

This range was found to vary with water depth, mixed layer depth, 
and sea surface temperature, The graphical representation can be seen 
in Figures 10 and 11 (note change in scale for range between Figure 
10 and Figure 11), These figures show the variation of range with sea 
surface temperature and mixed layer depth and with sea surface tempera- 
ture and water depth. 

Tie 100-decibel loss range can be determined using the following 


formula; 
D = 4000 
Bien = ie) - 0.12 mie (f= 10) + (2) -or (it) (5) 
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bP, 


(CB) 0525(T — 10) + aa (0,75) for M£& 100 meters, 
(II) = 0,155(T - 10) + aa eee (4,00) for M >100 meters; 
where 
R00 = distance to 100-decibel loss in nautical miles 
D = depth of water in meters 
T = sea surface temperature in degrees Centigrade 
M = mixed layer depth in meters, 
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CHAPTER V 


CONCLUSIONS 


iquation (4) was developed to predict the range to the starting 
point of the first three convergence zones and is accurate within 1500 
yards, Used in conjunction with Table 2 for convergence zone width and 
Figure 9 for convergence zone intensity, the operator can establish the 
location, width, and intensity of convergence zones one, two, and three, 
The field operator can obtain an accurate estimate of the 100-decibel 
propagation loss range by use of equation (5) where convergence zones 


forn, 
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CHAPTER VI 


RECOMMENDATIONS 


The convergence zone range should be examined more thoroughly 
for mixed layer depths between the surface and 100 meters, A more 
detailed examination over more shallow layer depths should permit 
refinement of the present formula for greater accuracy, 

Convergence zone width and convergence zone intensity should be 
investigated using a different method than the approach presented in 
this paper, A constant temperature profile might be maintained and 
the frequency of the source might be varied over a broad range of fre- 
quencies to find a relationship, 

The 100-decibel propagation loss range should be investigated 
using the excess depth below the convergence zone as a reference, 

More temperature profiles should be processed to refine the prediction 
formula at mixed layer depths other than 100 meters and for water 


depths other than 4000 meters, 
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APPENDIX 
Convergence zone range formula 
R= ([22,70+0,55(T-10)] Z 


+ [(r-.9) (0,2527+7?) - (2 - 1)%* (0,5)] BS ae 


This formula was developed in two parts, The first part related con- 
vergence zone range to sea surface temperature at the fixed mixed layer 
deptheef one hundred.meters, 

Using the least squares method a first-degree polynomial was 
fitted to the data for zone one, zone two, and zone three, The results 
of this fit gave the following polynomials; 

ZONE 1 22.70 wee 55¢T -"10) 

ZONE 2 45.399°+ 1 me2(T = 16) 

ZONE 3 68,07 + 1,644(T - 10) 


These three polynomials were then reduced to the polynomial: 


[22.70 + 0,55(T - 10)] z 


r 
I 


sea surface temperature in degrees Centigrade 


N 
i 


convergence zone number one, two, or three 


A plot of the actual data with the least squares curve fitted to 
it is given in Figure 12, As sea surface temperature increases it is 
necessary to consider deeper water depths in order for a convergence 
zone to forn, 

To develop the second part of the formula the first part was used 


as a standard, Then the change was calculated as a positive or 
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Figure 12 


“al: 


negative range as mixed layer depth increased or decreased from one 
hundred meters, This correction for the change in mixed layer depth 


is expressed by: 


ce Senne: 
M100 [ ( - 9) (0,25 24632) - (2 - 1)8 (0,5) 


M = mixed layer depth in meters 
(Note this entire term is equal to zero at M = 100), 
To represent the change in convergence zone range as a function 
of mixed layer depth and sea surface temperature, Figure 13 was developed, 
This figure shows the increasing divergence caused by changing mixed 
layer depth for convergence zones one, two, and three, This is a 
reduction of the data presented in Tables 4, 5, and 6, The following 


features should be noted from the graphical display: 
Jide — change in slope for convergence zones one, two, and three 


(Z - 1)9.8 (0,5) - term corrects for the initial range difference 
between the surface and 300 meter mixed layer 


depth 
 - term determines the sign of the correction, Thus an 


increase in range occurs for mixed layer depth greater 
than 100 meters and a decrease in range for mixed layer 
depth less than 100 meters, 


This formula is effective for sea surface temperatures between 
10 and 25C and for mixed layer depths from the surface to three 


hundred meters, 
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MLD = 000m 


actual range (nautical miles) 


predicted range (nautical miles) 


MLD = 100m 

P A 
22.70 23,00 
23,25 23800 
23,808 23,50 
Aigo >. 24800 
CaO 215 50 
2565 ie25)x50 
26,00 26,50 
20,9) eee 
27,10 27,00 
27.65 
28,20 28,00 
28.75 28,00 
2on00 29400 
29,85 29,50 
30,40 30,50 
30,92 

TABLE 4 


PREDICT#D AND ACTUAL RANGE OF CONVERGENCE ZONE ON: 


o_o See ews eee Sl ee 


MLD = 200m 


MLD = 300m 


22,86 
237 
24,30 
25,01 
296th 
26,45 
27 uae 
2 yc, 
28,60 
Agee 
30,04 
30 
31.16 
319 
32.90 
32 #@uL 





= predicted range (nautical miles) 


A = actual range (nautical miles) 


TCC) 


MLD = 000m 


Je 
45,36 
46.25 
A714 
48,03 
48,92 
49,82 
90,71 
51,60 
9249 
53.38 
54,18 
pe! 
56,06 
Jou? 
57 84 
28,73 


MLD = 100m 
P A 
45,40 46,00 
46,50 46,50 
L760, L760 
48.70 48ye@0e 
L780 = 49788 
5,92 50,58 
22 IO . 62,50 
23.40 |= FRO 
B20 53,50 
5530 
568h0 . BD, 50 
S120 56,50 
58,60 58,00 
5910 « So a8 
60,80 61,00 
61,90 
TABLE 5 
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PREDICTED AND ACTUAL 2ANGY OF CONVERGENCE ZONE TwO 


cn ee 


MLD = 200m 


MLD = 300m 


1 e, 
I 


> 
li 


mec ) 


See eC Se eS oe eee OES eee 


MLD = 000m 


68,04 
69.33 
70,63 
TADS 
ees 
Tee 
ore 
ieee. 
78,41 
79,71 
81,00 
82,31 
83,60 
84,90 
86,19 
87 49 


actual range (nautical miles) 


predicted range (nautical miles) 


MLD = 100m 
P A 
68,10 68,50 
69.75 70,00 
71 40 "71,00 
(Oe0>  “WeeO0 
Ge 70 "7350 
M6625 76,00 
TogOO I IOe 
(2g) cL. oe 
81,30 80,50 
82.95 
S960 "63850 
BOse> "4, 0U 
87.90 87,00 
ee Uo 
71,20 9 ilreo 
Saree 
TABLE 6 
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MLD = 200m 


MLD = 300m 


68,23 
(ORD8 
72.91 
Voea7 
77,64 
80,01 
SEO 
81, ,73 
87,08 
89 43 
91,80 
7 gee 
96,50 
98,85 


101.223 


113,57 
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